Multiple scaling regimes in simple aging models 
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We investigate aging in glassy systems based on a simple model, where a point in configuration 
space performs thermally activated jumps between the minima of a random energy landscape. The 
model allows us to show explicitly a subaging behavior and multiple scaling regimes for the corre- 
lation function. Both the exponents characterizing the scaling of the different relaxation times with 
the waiting time and those characterizing the asymptotic decay of the scaling functions are obtained 
analytically by invoking a 'partial equilibrium' concept. 
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The dynamics of glassy materials can be strongly de- 
pendent on the history of glass formation Gen- 
erally speaking, one finds that the relaxation dynamics 
becomes increasingly slower with the "age" of the sys- 

[ tem, that means with the time expired since the ma- 
terial was brought into the glassy state. Such aging phe- 
nomena have been identified in many systems and vari- 
ous dynamical probes (for a recent review, see e.g. [^). 
Prominent examples are shear stress relaxations in struc- 
tural glasses [Q, thermoremanent magnetizations in spin 
glasses 1^ , and electric field relaxations in dipolar glasses 
M. A convenient way to quantify aging in such exper- 
iments is to disturb the probe at time by a sudden 

' change of external field, and to measure the response 
R(tu, + t,tw) at a later time + t. Often, the charac- 
teristic relaxation time grows proportionally to the age 
t^. When t is larger than all microscopic times asso- 
ciated with fast time-translational invariant relaxations, 
one then expects R{tw + t,!^) to depend on the ratio 
t/t^ only, i.e. R{t^+t,t^) = F{t/t^). 

In principle, however, one cannot rule out other scal- 
ing forms, as e.g. R{tw+t,tw) = F{t/t!;^) with ^ > 

^ being different from one. In particular the case /i < 1, 
which has been called 'subaging' because the effective 
relaxation time grows more slowly than the age of the 
system, seems to be of experimental relevance jsj. It is 
moreover possible that there exist, for given waiting time 
iu,, various scaling regimes in time t, which are governed 
by different relaxation times oc s = 1,2,.... More 
precisely, depending on how t is scaled with t^, one can 
obtain different asymptotic scaling functions in the limit 
— > oo. For example, for /ii > /i2 > 0, one may find 

: i?(t»+AiC,^»)-^-F^i(Ai) andii'(t^+A2C,*»)~^2(A2) 
when — > oo. In fact, the occurrence of different scaling 
functions being associated with various time regimes has 
recently been conjectured on the basis of analytical re- 
sults for the Langevin dynamics of mean-field spin glass 
models |^J^. So far, however, it was not possible to val- 
idate these conjectures, or to exemplify them in some 
reasonable phenomenological models. 

In this Letter we will discuss a model that allows us 
to demonstrate for the first time explicitly the possi- 



ble occurrence of subaging behavior and multiple scaling 
regimes. This model, which has a strong resemblance to 
the previously studied "trap model" is motivated 

by the simple and widespread view that glassy dynam- 
ics may be described by a thermally activated motion of 
a point ("particle") that jumps among the deep (free) 
energy minima Ei of a complex configuration space. Ac- 
cording to extreme value statistics one may expect the 
distribution p{E) of these deep minima to be exponen- 
tial, and indeed, mean-field theories of spin glasses |ic| ] 
and recent results from molecular dynamics simulations 
[]ll| suggest this to be the case. 

To be specific, let us consider a d-dimensional cu- 
bic lattice and assign to each lattice site i an en- 
ergy Ei^ —oo < Ei < 0, drawn from the distribution 
p{E) = T~^ exp{E/Tg). The particle jumps among near- 
est neighbor sites only, and the jump rate from site i to 
a neighboring site j is 

w^j = lyexp {-l3[aEj - (1 - a)Ei]) , (1) 

where the "attempt frequency" iy = l sets our time unit, 
f3~^ = T is the temperature (or thermal energy), and 
the parameter a specifies how the energies of the ini- 
tial and target site are weighted. In order for the Wij 
to obey detailed balance, a can assume any real value, 
but on physical grounds it is reasonable to restrict a to 
the range < a < 1. Independent of a, the system un- 
dergoes a "dynamical phase transition" a,t T — T^: In 
the high-temperature phase, where T > Tg, the proba- 
bility (p{E) oc p{E) exp(— /3i?) for finding the system in a 
state with energy E is normalizable and thermal equilib- 
rium will be approached with a characteristic equilibra- 
tion time that diverges for T\Tg. By contrast, in the 
glassy phase, where T <Tg and ^p{E) is not normalizable, 
the dynamics never becomes stationary but 'ages'. 

It is important at this point to stress the differences 
between the above model and the earlier studied trap 
model In the latter, the jump rates depend only on 
the energy of the initial site corresponding to a = in 
eq. (|l|), and this allows a straightforward mapping onto a 
continuous time-random walk with a waiting time distri- 
bution decaying as a power law. Much more important. 
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FIG. 1. Double-logarithmic plot of Il{tw + t,tw), 1 — 
n(t™ + t,t^) (symbols o, □, a, V, o) and U{tu, + t,t^), 
1 ~'n{tw +t,tiu) (symbols >, +, x) as functions of t/t!^^ in 
d= 1 and d = 3 for (6*, a) — (1/4, 3/8). Different symbols cor- 
respond to different waiting times, t^ — W^ (o), 10^ (V), 10^ 
(a), 10» (□), 10^ (o) and U, = 10« (>), 10" (+), lO" (x) 
in d = 1, as well as U = lO" (o), 10^ (V), 10** (a), 10® (□), 
10"' (o) and U = 10" (>), lO" {+), 10^^ (x) in d = 3. The 
solid lines have slope e (small t behavior) and slope —S (large 
t behavior). For the value of the exponents e and 5, see text. 



the trap model was so far considered only on a mean- field 
level corresponding to an "annealed situation" , where the 
site energies are drawn anew after each jump. 

In order to study aging effects in the glassy phase, we 
focus, as in the trap model, on the (disorder averaged) 
probability Il(tyj+t,tw) that the system does not change 
its state between tyj and tw + 1. For particles hopping on 
a lattice, this can be interpreted as a dynamical struc- 
ture factor Initially, the particle is located on any of 
the sites and then it starts to explore the configuration 
space at some T<Tg. Physically, this means that we are 
considering an instantaneous quench from T = oo. 

Our idea to explore the scaling properties of 11 (i^ + 
t,tw) for both t and 1^ becoming large is based on 
the following "partial equilibrium" concept: From the 
time after the quench up to the particle has fol- 
lowed a Brownian path in configuration space that on 
average consists of S{tw) distinct and mutually con- 
nected sites. On a typical path with S ~ S{t^) sites 
then, it is reasonable to think that the particles should 
effectively equilibrate, i.e. the probability to be on a 
particular site j of the path may be approximated by 
Tj/Sfe=i ''"fci where Tj = exp(— /3i?j) (1 <Tj < oo). Con- 
ditioned on being at the site j, the system has a proba- 
bility exp{~tJ2n, Wj,n,) = e^Pi-tT^''^J2n, '^n,) tO 



change state within time t, where the sum over Uj runs 
over all nearest neighbor sites of j. Exactly two of these 
neighboring sites are considered to belong to the Brow- 
nian path in view of its one-dimensional topology. The 
remaining 2((i— 1) neighboring sites are assumed to have 
not been visited before. Hence, we may deduce the scal- 
ing properties of Il{tw+t, t^) from 



Tj exp -tT° 



(2) 



where (. . .) denotes an average over (2d — l)S{tm) uncor- 
related random numbers Tj that are distributed according 
to the power law (/)(r) = 0t~^~^ with e = T/Tg. 

Clearly, the partial equilibrium concept is an approxi- 
mation that needs to be tested. To this end we have de- 
termined Il{tw+t, tw) and S{tjju) in d= 1, 3 for various val- 
ues of 6 and a by Monte Carlo simulations. Then we took 
S{t^,) from these simulations to calculate Il{tm -\-t,tw) 
from eq. (^ . The disorder average in this simple numer- 
ical evaluation was performed over a set of realizations 
independent of the ones taken in the simulations. 

Figure |l| shows the results for one representative pa- 
rameter pair {9, a) — (1/4, 3/8). The data have been col- 
lected as functions of t for a broad range of fixed waiting 
times tw and are plotted already in scaled form as func- 
tions of Ai = t/t!^^ with exponents fj,i = fii{9,a) being 
specified below. As can be seen from the figure, the data 
for both Il{tw+t, tyj) and Il{t^+t, t^) collapse onto single 
curves i^i(Ai) and i^i(Ai), respectively. Although the 
two scaling functions are different, their asymptotic be- 
havior for large Ai is the same, Fi(Ai) ^ i^i(Ai) ~ Aj"* 
(see the solid lines in Fig. |]) . The values of the exponent 
5 = 6(6, a) are given below. 

It is important to inspect also more closely the 'short' 
time regime, where the complementary probability 1 — 
Ii{tw+t,tw) for the system to change state between 
and tw+t is small. This complementary probability can 
be as relevant as n(iu,-|-i, i^) dependent on the physical 
quantity being measured. Scaling plots of l—Ii{tyj+t, t^) 
and 1— n(t^,-|- 1, tyj) as functions of Ai = t/tj^^ in d= 1, 3 
are also shown in Fig. |l|. Again there is a good data col- 
lapse and for Ai ^ we find l-Fi(Ai) - l-Fi(Ai) - Af 
with exponents e = e(0,a) given below. An analogous 
overall behavior as displayed in Fig. ^ was found also for 
other pairs {9, a) (with < 6* < 1, < a < 1). Wc thus 
conclude that the partial equilibrium concept not only 
yields the correct scaling behavior (same /ii exponents) 
but also the correct asymptotics of the scaling functions 
|12|. However, a study of the 'participation ratios' (see 
1 13 1 for their definition) shows that the partial equilib- 
rium concept is not exact, even for large times Jl4]. 

Now we turn to the analytical study of eq. Since 
Il{tw+t, tw) depends on t^, only via S{tw) let us first dis- 
cuss the scaling of S(tw) with time tw For a = this 
problem has been addressed some time ago (see e.g. [p^) 
and for t ,„ — > oo one finds 
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(3) 



(In d = 2 there are logarithmic corrections, S{t^) ^ 
[tw / logt^]^ .) For < a < 1 one expects (||) not to change, 
since a only affects the nearest-neighbor hopping rates 
but not the transport properties on large length scales. 
In fact, in our simulations we always found eq. to hold 
true. Let us note also that in d = 1 one can give a sim- 
ple finite-size scaling argument to show that (j^) remains 
vahd for 0<a< 1. 

Next we derive the scaling properties of 11 as a function 
of t and S = S{tj^), and then use eq. (^) to obtain the 
corresponding scaling properties of 11 as a function of t 
and When replacing the denominator in eq. by 

dX exp(— A X)f=i '''k), the average over the tj can to a 
large extent be factorized, and one obtains the following 
asymptotic formula valid in the limit of large S 
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Here k = r(l— 6'), where r(.) is the Gamma function, and 
S = kS. Note that for t = and 5 ^ oo, eq. (§) yields 
the correct normalization IT— >1. 

After two transformations A SX^ and t S'^/^t we 
can identify Ai =t/5'(^~")/^ as a scaling variable corre- 
sponding to a first characteristic time ti'^ S(t^Y^~°'''/^ ~ 



t^l^ . We thus obtain 

7(1 - a) 

/^l = n 



(5) 



with 7 from eq. (|3|). This exponent /ii has been used to 
collapse the data in Fig. |l|, i.e. we took /ii — (1— a) / (14^) 
in d = 1 and /ii = (1 — a) in d = 3. Based on the equi- 
librium concept it is is easy to show that ti has a sim- 
ple physical interpretation: It scales as the typical max- 
imum trapping time tmax(itu) encountered after tw In 
d=l one thus finds a subaging behavior even for a = 
since the deepest trap is visited a large number of times 
N{t^)r^t'^J^^"^\ SO that ti-tmax(i»)'-^i»/A^(i»)<t«,. 
Similarly, for 0:7^0, the deepest trap is revisited a large 
number of times in all dimensions d due to a strong back- 
ward jump correlation when the particle leaves a site with 
low energy. 

The scaling function in the first time domain i = Ait^^ 
reads 
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FIG. 2. Double-logarithmic plot of 1 — Il{tw +t,tw) and 
1— n(fu,-|-t, tuj) as functions of t/t'^ in d=l and d = 3 for the 
same parameters as in Fig. ^ (the same symbols have been 
used for the various waiting times t^). The solid lines have 
slope one and 9 /a, and indicate the limiting behavior of F2(.) 
according to eq. (|ll|). 

and has the limiting behavior 
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The constants c>, c<, and can be expressed in terms 
of a, 6 and d but are not of interest here. Equation (|^) 
yields the exponents S = 9/{l — a) and e= (1 — 6')/(l — a) 
taken in Fig. 1 to characterize the decay of Ft and rise of 
1— n, respectively. For a = 0, one recovers the results of 
the annealed model namely S — 9 and e = l — d. 

For 6 > a the function Fi(Ai) describes the scaling 
properties of 11 completely. However, based on eq. (Q) 
one finds that for 9 < a < 1/2, there exists a second 
scaling variable A2 = t/^(i"2a)/e _ f/fX^i^^a)/^ ^ 
yielding 



7(1 - 2a) 



(8) 



Therefore, a second characteristic time scale t2 ^ 
5'(i-2a)/9^ diverging when tyj — > 00, governs the behav- 
ior when H is close to one. (Note that for fixed A2, 
Ai — 3^°"/^ K2 for S — > 00). In the new time domain 
one finds the following generalized scaling form. 



5(1-H) = F2(A2). 



'<a<l/2. 



(9) 
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where F2{.) is given by 



g{x) = xi- 



du 



u^+i 



— (1-e-^) exp[-(x/u)TT-^]. (10) 



Here -^oo =linir-»oo '4'{t)/{0t ^~^), and '0(t) denotes the 
probabihty distribution for the variable f = (r^+r^f )^/", 
i.e. V(T) = ((5(f-[rf +T2"]i/")). From (Il0|) foUows 
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where cq is a constant dependent on and a. We note 
that 1 — n matches continuously as one leaves the short 
time scaling regime (t ^ t!^^) described by F2 to enter 
the regime described by the scaling function Fi (where 
t ~ ) . Figure 2 shows the short time behavior of ft 
and n, rescaled as in (^ for the same parameters as in 
Fig. |l]. As can be seen from the figure, the data approach 
the two power laws predicted by eq. (|l^) for large tw 

For clarity, we illustrate the overall behavior of Il{tu]+ 
t,t^) as a function of t in Fig. 3. For {9, a) values in 
the two-timc-scaling region <9 <a<l/2 (shaded area 
of the a-0-diagram shown in the inset), there exist three 
different t regimes: (/) U{ti^+t,ti^) ^l—const. t~''[t/tl^^] 
for (//) n(t„ + i,t„)-l-const.'t-^[t/tJ^,^]«/" = 

1 - const.' [t/C]^/" for C <*«C> and (///) n{t^ + 
t,t^)r^[t/ti^^]-^ for t!^^ <i. When (61, a) lies in the un- 
shaded area of the a-6'-diagram, the first regime t ^ tl^^ 
does not exist (or, more precisely, it then becomes irrele- 
vant in the limit of large tyj). Note that for tw—^oo and 
^2 = t/t!ij^ fixed, the long-time regime in Fig. 3 "moves to- 
ward infinity" and the behavior is fully described by the 
second scaling function F2(A2), while for i^, — > 00 and 
Al = t/f^^ fixed the short-time regime in Fig. 3 "moves 
toward zero" , and the behavior is fully described by the 
first scaling function Fi(Ai). 

In summary we have shown (i) that generalized trap 
models can exhibit subaging behavior, induced by multi- 
ple visits to the same trap, and (ii) the possible existence 
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FIG. 3. Sketch of the behavior of n(i™-|-t, tm) as a function 
of time t in the three regimes {I-III). The shaded area in the 
Qf-S-diagram marks the two-time scaling region. 



of several distinct scaling regimes in the two-time plane. 
Such a possibility is of crucial importance for the inter- 
pretation of experiments, since the waiting time can typ- 
ically be varied between one minute and a few days only 
(with some notable exceptions [Q). The occurrence of 
several time regimes then may get masked by an appar- 
ent rescaling of the relaxation curves by a single effec- 
tive value of /z. From a theoretical perspective, it would 
be interesting to study the existence of a generalized 
Fluctuation-Dissipation Theorem in the aging regime, as 
it is predicted by mean-field spin-glass models ||l6| . This 
problem is intimately related to the validity of the partial 
equilibrium concept introduced above. 
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